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ABSTRACT 


To  implement  our  theoretical  method  for  defining  texture 
measurements*  it  ie  necessary  to  hove  the  abilityH*o  generate 

Z  „  * 

i  toituree.  In  this  report*  a  texture  generation^  method  is 

a 

* 

1  described  based  on  the  use  of  Markov  chains.  This  method  was 
implemented  using  a  Markov  Texture  Generation  Program. 
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1.  INTRODUCTION 


A  theoretical  .method  -  for  defining  texture  measurements  was 

ptc Vuu-S  v»rWJ  *  £, ! 

defined  in*-eefei  nrtH  1-  and"  2.  ■  To  implement  this  method*  it 

Cl 

is  necessary  to  have  the  ability  to  generate  textures.  Two 
texture  generation  procedures  are  proposed:  first*  a 

modified  Oagalowicz  generation  procedure  and  second*  a  Markov 
texture  generation  procedure^^Tl^.  The  objective  of  this 
report  is  to  describe  the  implementation  of  the  latter  of 
these  generation  procedures.  Before  describing  the 
implementation  of  the  Markov  tpxture  generator  *  it  is 
necessary  to  point  out  the  importance  of  the  texture 
generator  procedure. 


Texture  analysis  is  a  very  important  part  of  image 
processing.  Textural  feature  extraction  (measurements)  is  a 
key  point  in  the  conception  of  image  classification  schemes. 
Comparative  studies  of  the  efficiency  and  accuracy  of  those 
techniques  have  been  done  to  identify  their  advantages  C3*43 
using  real  image  data.  This  is  reasonable  because  the 
ultimate  purpose  is  to  analyze  real  data.  Conclusions  of 
such  studies  might  sometimes  be  not  general i zable  because  of 
the  limited  member  and  variety  of  the  data  samples.  In  order 
to  proceed  with  comparative  studies  in  a  more  through  manner* 
one  can  supplement  the  real  data  by  generating  texture 
patterns  having  specified  statistics.  The  character istics  of 
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•ach  technique  can  be  investigated  to  provide  stronger 

conclusions.  In  C53*  Harkov  generated  textures  are 

considered  for  a  theoretical  comparison  of  ’lour  texture 

algorithms.  They  examined  the  amount  of  texture-context 

Cl 

information  contained  in  four  algorithms.  Texture  generation 
also  has  important  applications  such  as  the  generation  of 
test  samples  for  the  analysis  of  the  human  visual  recognition 
system  or  for  the  control  of  texture  classification  algorithm 
C63.  Julesz  171  demonstrated  that  computer  generated  images 
with  controlled  statistical  topological  or  heuristics 
properties  are  liable  to  reveal  some  basic  organization 
principles  of  information  processing  in  the  sensory  nervous 
system.  These  patterns  deprive  subjects  of  their  life-long 
learned  habits  of  recognition  and  make  them  rely  on  more 
primitive  mechanisms.  Julesz  used  Harkov  textures  in  his 
investigation  of  human  texture  perception. 

The  generation  of  textures*  in  general*  can  be  divided  into 
two  schemes:  structural  analysis  CB-123  and  statistical 

analysis  C13-283.  Harkov  generation  procedure  belongs  to  the 
statistical  analysis  scheme. 

A  Harkov  texture  image  can  be  generated  using  the  outcomes  of 
gray  tone  values  from  a  regular  Harkov  chain  by  arranging 
these  outcomes  in  a  sequential  manner  along  an  image  line. 
Zn  a  more  formal  way*  the  generation  procedure  is  one  which 
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the  random  field.  X(n.  m)  »  0  £  n  <«*)  0  ^  m  0»i *  generated  a 

row  at  a  timet  each  row  ie  generated  in  a  lef t-to-r ight 

direction!  and  the  grog  level  of  the  left  most  pofnt  of  each 

>  row  is  asigned  according  to  an  initial  distribution.  The 
'  'j 
~  rest  of  the  points  of  each  row  are  assigned  grag  levels 

according  to  the  one-step  transition  probabi 1 ities.  p(i/j)» 

of  a  Markov  chain.  The  major  advantages  of  this  method  are 

a  simple  generation  process  and  a  direct  specification  of 

spatial  co-occurrence  matrix  of  generated  texture. 

Markov  generation  procedure  is  used  in  conjuction  with  design 
objective  2  Csee  reference  13.  This  is  needed  because  the 
desires  to  create  a  set  of  measurements  which  can  match  the 
capabilities  of  spontaneous  human  texture  perception  that 
"resembles"  the  primitive  mechanisms  of  human  perception,  and 
the  modification  requires  to  make  the  Cagalowicz  procedure 
severelg  restricted  the  class  of  texture  that  can  be 
generated  using  this  method. 

In  section  2  of  this  paper,  a  theorg  and  the  nonmeclature  is 
presented.  Section  3  describes  the  implementation  procedure. 
Finallg.  section  4  presents  some  sample  runs. 


2.  THEORY 


A.  Ointral  Theory 

Definition:  Markov  chain  C293.  A  Markov  chain  ie  aii  sequence 
of  experiments  performed  on  a  system  8  with  the 
following  properties: 

1.  At  any  given  time  the  system  is  in  one 

of  the  states  E1<E2* . En  and  the 

outcome  of  each  experiment  is  that  8 
either  is  unchanged  or  is  changed  to  a 
new  state  from  the  set  <E1<E2> . En>. 

2.  The  system  can  change  from  one  state  Ei 
to  another  state  Ej  only  as  a  result  of 
the  experiment. 

3.  The  probability  of  changing  the  system 
from  the  state  Ei  to  the  state  Ej 
depends  only  on  Ej  and  on  Ei«  and  the  state 
of  the  system  at  the  time  of  the 
experiment.  This  probability  is  denoted 

by  p(j/i>. 


In  a  more  formal  way «  consider  a  stochastic  process 
X«<Xn;neN>  with  a  countable  state  space  E.  A  Markov  chain  is 
a  sequence  of  random  variables  such  that 

foU  «  E 
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That  i«>  the  next  state  Xn+1  is  independent  of  the  pest 
states  XO. . Xn-1  provided  that  the  present  state  Xn  can  be 

known.  Since  p(j/i)  gives  the  probability  for  thf  transition 

r  ‘ 

r  of  S  fro*  state  Ei  to  the  state  Ej#  it  is  called  the 
^  transition  probability  for  the  Markov  chain.  It  is  customary 
to  arrange  the  p<j/i)  into  a  square  array  and  to  call  the 
resulting  matrix  Cp<j/i)3  the  transition  matrix  of  the  Markov 
chain.  The  transition  matrix  has  the  following  properties: 

1.  The  transition  matrix  represents  probabilities. 
The  elements  of  this  matrix  p(i#j>  are  all 
nonnegative. 

£(.i/i)>,o  -fan  £  E. 

2.  For  each  row  of  the  matrix  the  sum  of  the 
elements  in  that  row  must  be  1. 

mo  i  { $  e. 

Any  matrix  that  satisfies  these  two  condition  is  called  a 
stochastic  matrix. 

Theorem  C303:  For  any  n»  m#  e  N  with  m  )  1  and  iO,  .  ...  im  6  E. 

. 

Corollary  £303:  Let  TTo  be  a  probability  distribution  on  E» 

and  suppose  p<Xo«i>"TTo( i )  for  all  i  (  E. 

Then  for  any  m  f  N  and  io* . . .  .  im  #  E , 

£  U-i/ta) . -l) 

This  corollary  shows  that  the  joint  distribution  of 
Xo»Xl»....Xm  is  completely  specified  for  every  m  once  the 
initial  distribution  TTo  and  the  transition  matrix  CP3  are 


known. 


B.  Method  of  Generating  Texture 


The  textures  to  be  generated  are  ones  whirfi  can  be 
represented  by  random  fields  generated  using  Markov  chains. 
It  generates  a  row  at  a  time)  each  row  is  generated  in  a 
lef t-to-r ight  directions!  and  the  grey  level  of  the  left-most 
point  of  each  row  is  assigned  according  to  an  initial 
distribution  TTo.  The  rest  of  the  points  of  each  row  are 
assigned  gray  levels  according  to  the  one-step  transition 
probabilities.  p(j/i).  of  a  Markov  chain.  Using  this 


ition  procedures  the  only  parameters  neccesary  to 


completely  specify  the  random  field  X(n>m)  are  the  initial 
distribution  TTo(i>.  O  .£  i  £  1-1.  and  the  one-step  transition 
probabilites  p(j/i>#  0  .£  i  0  <  4  £  1-1  ..where  1  is 

the  number  of  gray  levels.  To  facilitate  the  computation  of 
the  required  expected  values.  only  Markov  chains  which 
satisfy  the  following  three  constraint  conditions  are  used:  \ 

1.  Markov  chains  which  have  stationary  one-step 
transition  probabilities 

4>  (jit)  ,  O&Cii-l  ,  0< 

2.  Markov  chains  which  have  an  uniform  stationary 

distribution  TTs.  i.  e. . 

,  oi  tti-/ 

TT$  (oj  «  O  - 

3.  The  initial  distribution  must  be  the  uniform 
distribution,  i. e. » 

TToLO*  ili  ,  i*A-i 
V.  (0  ^  O  • 
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The  third  condition  guarantors  that  TTo  *  TTs.  Tho  resulting 
two-dimensional  random  fields  have  the  following  properties:  j 

1.  For  any  n#m#k  and  1>.  n*k.  the  random  Variables 

r  ‘ 

* 

X(n>m)  and  X(k<l>  are  statistically  independent. 

Ci 

2.  For  any  n»  m,  and  1«  the  random  variables  X(n>m)  j 

and  X(k«l)  are  related  by  the  m-1  step  transition 
probabilities.  These  transition  probabilities  can 

be. computed  directly  from  the  one-step  transition 
probabilities  of  the  Markov  chain  used  to  create 
the  random  field. 

3.  The  random  field  X(n»m)  is  stationary  with 
respect  to  arbitrary  trans lations. 


The  following  equation  relates  the  Markov  chain  parameters 
used  to  generate  a  texture  with  the  spatial  gray  level 
dependence  matrices  which  can  be  computed  from  this  texture. 

EjsU,i|S)j  *  |  *'z/- ••• 

[TToU)  TTo(<i>)  ult SAilJyJAt. 

where  p  <j/i>  is  the  d-step  transition  probability  of  the 


Markov  chain. 


3.  IMPLEMENTATION 


The  subroutine  MRKV  roads  from  a  given  data  file  the  number 
of  elements  in  state  space  and  the  block  size  of  each 

pattern.  The  texture  options  (see  fig.  i>  are*-  promted  to 

r  * 

the  CRT  to  be  entered  by  the  user.  The  state  sp4£e  patterns, 
transition  matrices*  and  initial  distribution  matrices  are 
read  from  the  given  data  file  in  that  order.  The  number  of 
gray  levels  are  equal  to  the  number  of  elements  in  state 
space.  Using  the  number  of  gray  levels  and  the  gray  level 
options*  the  program  computes  the  shade  values.  The 
subroutine  GMRKV  is  called  to  generates  the  texture. 

Subroutine  GMRKV  generates  a  rout  in  a  lef  t-to-r  ight 
direction.  The  grey  level  of  the  left-most  point  of  each  row 
is  assigned  according  to  an  initial  distribution  and  a 
uniform  random  number  generator.  The  rest  of  the  points  of 
each  row  arc  assigned  grey  level  according  to  the  transition 
probabilities  (TRANS)  and  an  uniform  random  number  generator 


Others  auxiliary  subroutines  are  needed.  Subroutine  RANDOM 
generates  a  real  random  number  between  0  and  1.  0.  Subroutine 
RPALOC  and  CPALOC  are  used  to  allocate  an  output  file  in  ELAS 
format. 


B.  User's  Guide 


To  generate  a  Markov  texture,  the  user  has  to  be  loaded  in 
ELAS.  A  control  file  needs  to  be  created.  The  user  has  to 
assign  ID1  to  the  input  file  that  contains  the  information 

X  t 

r 

necessary  to  generate  the  texture.  The  overlay  MKRV  is 

‘4 

yanked  and  the  procedure  necessary  to  generate  the  texture  is 
started.  An  example  is  presented  as  follow: 

♦ELAS 

-what  control  file? 

->MRKV.  CF 
-file  mgr? 

— >AC  MARKOV.  DAT  ID1 
-file  mgr? 

->LF 

1  idl  markov.  dat 
-file  mgr? 

— >YANK  MRKV 

-choose  one  of  the  following  texture  options: 

-  id«0  whole  image  same  texture 

-  id*l  left  half  one  texture* right  half  another  texture 

-  id-2  two  textures*  one  in  a  small  box  on  image 

-  id *3  top  half  one  texture*  bottom  half  another  texture 

->1 

-choose  the  gray  level  options: 

-  ig-0  assign  0  as  black 

-  ig*l  assign  0  as  white 


>0 


-image  file  needs  to  be  allocated 

—input  name  of  file  to  be  allocated 
->HARKOV.  IMG 

—read  the  block  size  (n  x  n)  of  each  pixel 

— >2 

-markov  generation  program  termination 
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-file  mgr? 

->LF 

1  idl  markov.  dat 

2  mktx  markov.  img 
-file  mgr? 

To  display  the  texture  generated*  follow  the  instructions  to 
display  an  image  given  in  ELAS  manual.  In  the  example  above 
all  the  capital  letter  instructions  are  given  by  the  user. 
The  lower  case  letter  instructions  are  promted  by  the  program 
on  the  CRT. 


The  input  file  containing  the  information  necessary  is  create 
using  the  Editor.  The  format  is  as  follow: 

1.  Number  of  elements  state  space  and  block  size 
of  each  patterns.  Interger 

2.  State  space  patterns.  Interger 


3.  Transition  matrices.  Real 


4.  Initial  distribution.  Raal 


An  aiampla  of  bow  the  data  file  has  to  be  created'^to  generate 


r  i 


a  Markov  texture  is  given  as  follow: 

2  2  number  of  elements  state  space  and  bloejt  size 


» »  1 

22  I  state  space 

2  2  J 
0.0  1.0  ^ 

1.0  0.0  l 

O.  5  0.  5  ( 

O.  5  0.  5  J 


patterns 


0.33  0.33 
0.33  0.33 


} 


transition  matrices 

initial  distribution 


Appendix  A  presents  a  listing  of  the  program.  Next  section 
presents  several  texture  generates  using  this  program. 
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RESULTS 

T«bl«  1  shows  ths  transition  matrices  anil  :th«  initial 

*  r  j 

distribution  matrices  for  texture  A  and  B.  Figuras  2  show  a 

-*  * 
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texture  generated  using  the  transition  matrix  for  texture  A * 
with  three  space  state  patterns  of  a  block  size  of  two*  three 
gray  level  (O  is  black)*  and  a  pixel  size  of  four. 


(A) 


(B) 


.» 

0 

.2 

.2 

.6 

.2 


0 

.4 

.6 

.6 

.4 

0 


.2 

.6 

.2 

.2 

0 

.1 


(C)  11/3 


1/3  1/31 


Table  1.  The  Markov  chain  parameters  used  to  generate 
the  textures.  (A)  The  one-step  transition  ma¬ 
trix  used  to  create  the  texture  on  the  left. 
(B)  The  one-step  transition  matrix  used  to 
create  the  texture  on  the  right.  (C)  The  sta¬ 
tionary  distribution  of  both  Markov  chains. 
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Figurs  3* 4# 5* 6  show  texturs  pairs  generated  using  the 

transition  matrices  for  texture  A  and  B  with  three  space 

state  patterns  of  a  block  size  of  two*  three  gray  level  (0  is 

black)*  and  a  pixel  size  of  one*  two*  three*  and  four 

respectively.  The  texture  option  was  one. 


fol 


Figure  2.  Texture  generated  using  the  one-step  transition 
aiatrix  (A) ,  and  the  stationary  distribution  ma¬ 
trix  (C) . 
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Figure  3.  Texture  pair  generated  using  the  one-step  tran 
sition  matrices  (A)  and  (B),  and  the  stationary 
distribution  matrix  (C) .  Pixel  site  is  one. 


■  •  •  *  j»  '  jk  .  *•» 

.'A’-.V-.'A'  k'  WV.  •  •, 


"7Tn»r*  YwW  V  \*  “j* 


Figure  5.  Texture  pair  generated  using  the  one-step  tran¬ 
sition  matrices  (A)  and  (B) ,  and  the  stationary 
distribution  oat r lx  (C).  Pixel  size  is  three. 


Figure  4.  Texture  pair  generated  using  the  one- step  tran¬ 
sition  matrices  (A)  and  (B) ,  and  the  stationary 
distribution  matrix  (C).  Pixel  size  is  two. 
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Figure  6.  Texture  pair  generated  using  the  one- step  tran¬ 
sition  matrices  (A)  and  (B) ,  and  the  stationary 
distribution  matrix  (C) .  Pixel  size  is  four . 


Figure  7  end  8  show  texture  pairs  with  the  same 
charac taristics  of  the  other  texture  pairs  except  that  the 
pixel  size  is  two  and  the  texture  options  are  two  and  three 
respectively.  All  this  figures  demonstrate  the  ability  of 
our  Markov  Teture  Generation  Program. 
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Figure  7.  Texture  pair  generated  using  the  one-step  tran¬ 
sition  matrices  (A)  and  (B) ,  and  the  stationary 
distribution  matrix  (C) .  Pixel  size  is  two.  Tex¬ 
ture  option  is  2. 


Figure  8.  Texture  pair  generated  using  the  one-step  tran¬ 
sition  matrices  (A)  and  (B) ,  and  the  stationary 
distribution  matrix  (C).  Pixel  size  is  two.  Tex¬ 
ture  option  is  3. 
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TMDS  SUBROUTINE  MRKV 


Language: 


Purpose  : 


In/Output 


Fortran  VII*  using  TMDS  subroutines;  for 
Perk  in-Elmer  8/32. 


To  handle  the  input  data  to  obtain  the 
information  necessary  to  generate  Markov 
Textures. 


Respectively*  the  number  of  elements 
state  space*  block  size  of  each  patterns 
state  space  patterns*  transition 
matrices*  and  initial  distribution. 


Usaoe 


CALL  MRKV 


uuuuou  uuu  uuu  uuuo 


FILENAME  -  LSU1 : MRKV.  FTN 


1 

2 

3 

4 
9 
6 

7 

8 
'9 
10 
11 
12 

13 

14 
19 
16 

17 

18 

19 

20 
21 
22 

23 

24 
29 
26 

27 

28 
29 


MATCH 

C 

C 

C  TEXTURE  MEASUREMENT  DEFINITION  SYSTEM 

c  i 

X  REMOTE  SENSING  AND  IMAGE  PROCESSING  LABORATORY  1 

t  : 

X  ELECTRICAL  AND  COMPUTER  ENGINEERING  DEPARTMENT 

G  LOUISIANA  STATE  UNIVERSITY.  BATON  ROUGE.  70883 
C 

C - 

C 

C  PROGRAM  MRKV  . 

C 

C 

C  AUTHOR 
C 

C  R.  E.  VASQUEZ-ESPINOSA 

C 

C  DATE 
C 

C  OCTOBER  1982 

C 

C - 

SUBROUTINE  BEOIND 
CALL  MRKV 
RETURN 
END 


SUBROUTINE  MRKV 


SUBROUTINE  MRKV 
INTEGER  STATES ( 32*  32/ 8).  CB(  10) 

DIMENSION  TRANS(8, 8. 2),PR0B(8,  2>» ISHADE<64) 

DATA  INPUT/ 'ID1'/,NCFSF/'FMQR'/ 

GET  LOGIC  UNIT  FOR  THE  INPUT  FILE 

CB( 1 ) -INPUT 
CALL  CFSUB(8/CB/CB) 

IF<CB<2)  .  LE.  0)  GO  TO  2000 
LUS-CB<2> 

GET  LOOIC  UNIT  FOR  THE  PRINT  FILE 
CALL  PRF(LUPRF) 

READ  NUMBER  OF  ELEMENTS  IN  STATE  SPACE  (NSTATE)  AND 
READ  BLOCK  SIZE  OF  EACH  PATTERN  (NSIZE) 

99  READtLUS,  * INSTATE.  NSIZE 

96  WRITEILUPRF,  100)  NSTATE.  NSIZE 

97  100  FORMATdOX.  'NSTATE-  I4.2X.  'NSIZE- 


14,  /) 


KRnranmiuv1  pirorrr  cro  •.*  v-  iw  *.-  % 


C 

C 

c 


c 

:c 

'c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


CHOOSE  TEXTURE  OPTIONS 

WRITE(6.  101)  f 

READ  (5,102)  ID  I. 

DETERMINE  NUMBER  OF  TEXTURES  INVOLVED  - 

IF(  ID  .EQ.  0)  NTEXT-1 
IF(  ID  .  QT.  0)  NTEXT-2 

READ  STATE  SPACE  PATTERNS  (STATES)  FROM  THE  INPUT  DATA. 

DO  10  L-l,  NSTATE 
DO  10  1*1,  NSIZE 

READ ( LUS,  * )  (STATES (I,  J,  L>  ,  J-l.  NSIZE) 

WRITE(LUPRF,  103)  (I,  J,  L,  STATES (  I,  J,  L) ,  J-l,  NSIZE) 

103  FORMAT ( 10X,  'STATES  ( '.  12,  ',  ',  12,  ',  12.  ')-  ',  14) 

10  CONTINUE 

CHECK  FOR  ZERO  CONTENT. 

DO  20  L-l, NSTATE 
DO  20  1*1,  NSIZE 
DO  15  J*l. NSIZE 

IF(STATES(  I,  J,  L) )  15,30,15 
15  CONTINUE 

20  CONTINUE 
GO  TO  40 

30  DO  50  L*l,  NSTATE 

DO  50  1*1, NSIZE 
DO  50  J-l,  NSIZE 

STATESd,  J,L ) -STATES (I,J,L)  +  1 

50  CONTINUE 
40  CONTINUE 

READ  TRANSITION  MATRICES  (TRANS)  FROM  THE  INPUT  DATA  FILE. 

DO  60  L-l.NTEXT 

DO  60  1*1.  NSTATE 

READ  (LUS,  *)  (TRANS  ( I.  J,  L),  J-l.  NSTATE) 

WR I TE (LUPRF,  104 )  (  I ,  J,  L,  TRANS ( I,  J,  L),  J-l,  NSTATE) 

104  FORMATdOX,  'TRANS (  12,  ',  12,  '.  ',  12,  ')-  ',F10.  5) 

60  CONTINUE 

DO  70  L-l.NTEXT 

DO  70  1*1 , NSTATE 
DO  80  J*2» NSTATE 

TRANSd,  J,L  > -TRANS  (I.J.L)  ♦  TRANS(  I,  J-l,  L> 

80  CONTINUE 

TRANiBd,  NSTATE,  L)-l.  0 

WRITE(LUPRF,  104)  (I.  J,  L,  TRANSd,  J,  L),  J-l,  NSTATE) 

70  CONTINUE 

READ  THE  INITIAL  DISTRIBUTION  (PROB)  FROM  THE  INPUT  DATA  FILE. 
DO  90  L-l.NTEXT 

READ  (LUS,  *>  (PROB(  J,  L),  J-l,  NSTATE) 


UUP  uuu  uuu  uuu  uuu  uuu 


WRITE(LUPRF» 105) !  J,  L,  PROB ( J,  L ) ,  J-l ,  NSTATE ) 

105  F0RMATU0X,  'PROB!  M2,  '.  M2.  ')«  '.F10.5) 

90  CONTINUE 

DO  95  L*l»  NTEXT 

DO  94  J-2, NSTATE  - 

PR0B!J,  L)-PR0B!U,  L)  +  PROB! J-l,  L> 

94  CONTINUE 

WRITE!LUPRF»  105)  !J,  L,  PROB!  J,  L)»  J«l»  NSTATE) 

95  CONTINUE 

CHOOSE  THE  CRAY  LEVEL  OPTION  AND  COMPUTE  THE  SHADE  VALUES. 

NGRAY-NSTATE 
WRITE (6/  106) 

READ  (5* 107)  10 

INITIAL  VALUE  FOR  SHADE  VALUES. 

N-NORAY  -1 

I SHADE ( 1 )— 0 

IF( 10  .  EQ.  0)  IC-1 

IF(  10  .EQ.  1)  IC— 1 

IF<  10  .EQ.  1)  I  SHADE  ( 1 )  — N 

DO  96  I— 1,  N 

I SHADE  (I  +  1)-  ISHADE  (I)  +  IC 
112  FORMAT!  10X.  'ISHADE!  '»  13.  ')«  M4> 

96  CONTINUE 

COMPUTE  SHADE  VALUES. 

FACTOR-255.  0/N 
DO  97  1-1 .  NORAY 

ISHADE! I ) -ISHADE  !I)  *  FACTOR 
WRITE!LUPRF»  112)1.  ISHADE!!) 

97  CONTINUE 

OENERATE  IMAGE. 

CALL  GMRKV!NSTATE,  NSIZE, PROB. TRANS, ID, STATES, ISHADE) 

WRITE!6.  108) 

RETURN  CONTROL  TO  FILE  MANAGER  !FMOR). 

CB! I )-NCFSF 

CALL  CFSUB ! 1 1 ,  CB,  CB ) 

RETURN 

2000  WRITE!6,  110) 

•FORMAT  AND  ERROR  STATEMENT. 

101  FORMAT! 'I.CHOOSE  ONE  OF  THE  FOLLOWING  TEXTURE  OPTIONS:  / 

*'li  ID-  0  WHOLE  IMAGE  SAME  TEXTURE',/ 

*'6  ID-  1  LEFT  HALF  ONE  TEXTURE,  RIGHT  HALF  ANOTHER  TEXTURE'./ 

•  'Ii  ID-  2  TWO  TEXTURES,  ONE  IN  A  SMALL  BOX  ON  IMAGE'./ 

*'li  ID-  3  TOP  HALF  ONE  TEXTURE.  BOTTOM  HALF  ANOTHER  TEXTURE',/ 

•  ) 

102  FORMAT! ID 


172  106  FORMAT!  '&CHOOSE  THE  CRAY  LEVEL  OPTIONS:  ',  / 

173  *'».  IQ«  0  ASSION  0  AS  BLACK'*  / 

174  #'*  10-  1  ASSION  0  AS  WHITE'*/) 

175  107  FORMAT  (ID  ? 

176  108  FORMAT <  '&MARKOV  GENERATION  PROORAM  TERMINATION.'/) 

177  ^  110  FORMAT  (  'fc****ERROR  WITH  LOOIC  UNIT  FOR  INPUT  FILE****',/) 

178  5  RETURN  C. 

179  ’S  END 

180  *»BEND 


v- .  •  V*  .v.^ v,v.-w  y."  v 


TMDS  SUBROUTINE  GMRKV 


Language  Fortran  VII#  using  TMDS  subroutine!  for 

Perkin-Elemer  8/32. 


Purpose  To  gonorato  Markov  textures. 


In/Output  :  Respectively#  the  input  from  the  subroutine 

MRKV.  The  output  is  a  contiguous  file  that 
contains  the  texture  generated  in  Elas 
format. 


Usage  CALL  QMRKV<NST ATE.  NS I ZE,  PROB.  TRANS.  ID.  STATES. 

IDSHADE). 

NSTATE  -  input  state  space  elements 
NSIZE  -  input  block  size  of  patterns 
PROB  -  input  initial  distribution. 

TRANS  -  input  transition  matrices. 

ID  -  input  texture  options. 

STATES  -  input  state  space  patterns. 

IHSHADE-  input  gray  level  values. 


Program  Logic  : It  generates  a  row  at  a  timei  each  row  is 
generated  in  a  lef t-to-r ight  directions! 
and  the  grey  level  of  the  left-most  point 
of  each  row  is  assigned  according  to  an 
initial  distribution.  The  rest  of  the 
points  of  each  row  are  assigned  gray 
level  values  according  to  the  one-step 
transition  probabilities  (TRANS>. 


FILENAME  »  LSU1: CMRKV.  FTN 


1  C 

2  C 

3  C  TEXTURE  MEASUREMENT  DEFINITION  SYSTEM 

4  C 

9  C  REMOTE  SENSING  AND  IMAGE  PROCESSING  LABORATORY  ^ 

6  3c 

7  IC  ELECTRICAL  AND  COMPUTER  ENGINEERING  DEPARTMENT  ** 

8  re 

*9  'C  LOUISIANA  STATE  UNIVERSITY.  BATON  ROUGE.  70803 

10  c 

11  c - 

12  C 

13  C  PROGRAM  OMRKV 

14  C 

15  C 

16  C  AUTHOR 

17  C 

18  C  R.  E.  VASGUEZ-ESP I NOS A 

19  C 

20  C  DATE 

21  C 

22  C  OCTOBER  1982 

23  C 

24  C - 

29  C 

26  C 

27  C 

28  SUBROUTINE  OMRKV ( NST ATE, NSIZE,  PR OB,  TRANS,  ID,  STATES, 

29  * I SHADE) 

30  INTEGER  STATES (32,  32,  8) 

31  DIMENSION  IB(912),  TRANS (8,  8,  2>,  PR0B(8, 2>,  KBUFF(512> 

32  * I SHADE (64) 

33  L00ICAL*1  PIXEL (912),  LBYT(4) 

34  EQUIVALENCE  (IP.LBYT) 

39  DATA  NAME/'MKTXV 

36  C 

37  C 

38  CALL  PRF(LUPRF) 

39  C 

40  C 

41  WRITE(LUPRF, 9300)NSTATE,  NSIZE 

42  9300  FORMATdOX.  'NSTATE-  M3,  2X,  'NSIZE-  M3) 

43  DO  9100  L-l, NSTATE 

44  DO  9100  1*1, NSIZE 

49  DO  9100  J-l,  NSIZE 

46  WRITE(LUPRF,  9200)  I ,  J,  L,  STATES ( I , U, L ) 

47  9200  FORMATdOX,  'STATES!  ',  12,  ',  ',  12,  ',  12,  ')  *'.  14) 

48  9100  CONTINUE 

49  C  INITIAL  VARIABLES 

50  c 

81  IL*  1 

92  LL*  912 

S3  IE*  1 

94  L£-  912 

99  MBR*  912 

NC*  1 

87  IRAN-  123496789 


98 

C 

99 

C 

CHECK  IF  THE  IMAGE  FILE  IS  ALREADY  ALLOCATED. 

60 

c 

6t 

IB(1)-NAME 

62 

CALL  CFSUBC8* IB* IB) 

63 

*w 

LU-IB(2) 

64 

IF(LU  .OT.  0)G0  TO  200 

69 

■*3» 

WRITE(6.  100) 

~66 

* 

IB(1)«NAME 

67 

- 

IB(3)-IL 

68 

1B(4)-LL 

69 

IB(9)-IE 

70 

IB(6)-LE 

71 

IB(7)«NC 

r 

72 

CALL  RPALOC(LU,  NAME. IB) 

73 

IF(LU  .LT.  0)00  TO  9000 

74 

IB(1)»  NAME 

79 

IB(4)«  LL 

76 

IB(3)«  IL 

77 

IB(9)«  IE 

78 

IB(6)«  LE 

79 

IB(7>—  NC 

80 

CALL  ELTRAN  (LU,  4*  0.  200.  IB* 1ST) 

81 

GO  TO  400 

82 

c 

*83 

c 

FILE  WAS  PREVIOUSLY  ALLOCATED*  CHECK  FOR  CORRECT  J 

84 

c 

89 

200 

CALL  ELTRAN  (LU.  3.  0.  200.  IB.  1ST) 

86 

IF( IB(1 )  .  NE.  NAME)  GO  TO  9002 

87 

IF(  IB(3)  .NE.  IL)  GO  TO  9004 

88 

IF( IB(4)  .NE.  LL)  GO  TO  9006 

89 

IF( IB(9)  .NE.  IE)  00  TO  9008 

90 

IF(  IB(6)  .NE.  LE)  GO  TO  9010 

91 

IF( IB(7)  .NE.  NC)  00  TO  9012 

92 

4 00 

CONTINUE 

93 

c 

94 

c 

LET  READ  THE  BLOCK  SIZE  OF  EACH  PIXEL 

99 

c 

96 

WRITE (A,  101) 

97 

READ  (9. 102)  NN 

98 

11*0 

99 

IN-912/NN 

100 

DO  30  IROW-1,  IN 

101 

* 

IF(I1  .  EQ.  NSIZE)  1 1«0 

102 

11-11  ♦  1 

103 

IF( 1 1  .NE.  1)  GO  TO  40 

104 

c 

109 

c 

DETERMINE  CURRENT  STATE  FOR  IMAGE  GENERATION. 

106 

c 

107 

IX- IRAN 

108 

CALL  RANDOM  ( IX. IRAN.  XRAN) 

109 

00  90  I— l.  N8TATE 

110 

IF  (XRAN  .LE.  PROB(I.D)  GO  TO  60 

111 

90 

CONTINUE 

112 

60 

ICURNT-  I 

112 

NS-  NSIZE  •  NN 

114 

DO  70  J-l.  912.  NB 

M 


119 

JCOL-  J 

116 

C 

117 

C 

DETERMINE  TEXTURE  INDEX.  + 

118 

c 

119 

INDX-  1  ^ 

120 

IF  < ID  .  EQ.  0)  00  TO  71 

121 

IA  -  IROW  *  NN 

122 

** 

IF(  ID  .EQ.  1  .AND.  JCOL  .  OT.  256  )  00  TO  72 

123 

.  . 

IFCID  .EQ.  3  .AND.  IA  .  OT.  256)  GO  TO  72 

124 

* 

IF( ID  .  NE.  2)  00  TO  71 

129 

IF(  ID  .  LE.  320  .OR.  IA  .  OT.  448)  GO  TO  71 

126 

IF(  JCOL  .LE.  320  .OR.  JCOL  .  OT.  448)  00  TO  71 

127 

72 

INDX-2 

128 

71 

CONTINUE 

129 

c 

130 

c 

DETERMINE  CURRENT  STATE. 

131 

c 

132 

IX  -  IRAN 

133 

CALL  RANDOM  <  IX,  IRAN.  XRAN) 

134 

DO  80  K-  1 , NSTATE 

139 

IF  <  XRAN  .LE.  TRANSdCURNT,  K,  INDX))  00  TO  90 

136 

80 

CONTINUE 

137 

90 

ICURNT-  K 

138 

KBUFF(J)-  K 

139 

70 

CONTINUE 

140 

c 

141 

c 

OEN 

"RATE  EACH  LINE  (912  PIXELS). 

142 

c 

143 

40 

DO  209  JLK-1.912.NB 

144 

JMAX-  JLK  ♦  NB  -  1 

149 

IF( JMAX  .OT.  912)  JMAX-912 

146 

J2-  NN 

147 

Kl-  KBUFF(JLK) 

148 

DO  210  K-  JLK,  JMAX 

149 

J1-J2/NN 

190 

CCCCCCCCCCCCCCCWR X TE CLUPRF , 9200)11.  Jl,  Kl, STATES ( 11. Jl, K1 >CCCCCCCCCCCCCCC 

191 

M-  8TATES  ( 11.  Jl. Kl ) 

192 

IP-  I8HADE  (M) 

193 

CCCCCCCCCCCCCCCWR ITC(LUPRF< 111 )CIR0W.  JLK,  K, M, IPCCCCCCCCCCCCCCCCCCCCCCCCC 

194 

CC1UCCCCCCCCCCF0RMAT(9X.  'IROW-C  '.  13.  2X.  'JLK-C  13.  2X.  'K-C ',  13,  CCCCCCCCC 

199 

CCCCC*CCCCCCCCC2X,  'M*C  * •  13.  2X,  'IP-C I4)CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

196 

PIXEL(K)-  LBYT (4 ) 

197 

J2-  J2  ♦  1 

198 

210 

CONTINUE 

199 

209 

CONTINUE 

160 

c 

. 

161 

C 

PUT 

THE  LINE  IN  THE  OUTPUT  FILE. 

162 

c 

163 

NNX-NN 

164 

DO  220  J  -  1,  NN 

169 

1  *  IROW  •  NN  -  NNX  ♦  1 

166 

CALL  RDWR(LU.  4,  1, 1,  IB,  PIXEL,  1) 

167 

NNX*  NNX  -  1 

168 

220 

CONTINUE 

169 

30 

CONTINUE 

170 

RETURN 

171 

c 

BMataa&aaaa^ 


■ .  i  , 


^’1 


UU  i . •  CMJ  O 


ERROR  STATEMENTS. 

9000  WRITE(6,  9001 ) 

9001  FORMAT  (  '«c**#ERR0R  IN  ALLOCATING  THE  OUTPUT  FILE'#/) 
RETURN 

9002  WRITE(6, 9003) 

9003  FORMAT  <  'lc###ERROR  USAGE  NAME  NOT  MATCH',/)  - 
RETURN 

9004  WRIT£(6,  9005) 

9003  FORMAT  (  'Sc**#ERROR  INITIAL  LINE  NOT  MATCH'./) 

RETURN 

9006  WRITE (6.  9007) 

9007  FORMAT < '&*##ERROR  LAST  LINE  NOT  MATCH',/) 

RETURN 

9008  WRITE (6,  9009) 

9009  FORMAT < 'fc*««ERROR  INITIAL  ELEMENT  NOT  MATCH',/) 

RETURN 

9010  WRITE(6,  9011 ) 

9011  FORMAT  <  'Jt***ERROR  LAST  ELEMENT  NOT  MATCH',/) 

RETURN 

9012  WRITE(6#  9013) 

9013  FORMAT  <  'lc***ERROR  CHANNEL  NUMBER  NOT  MATCH',/) 

RETURN 

FORMAT  STATEMENT. 

100  FORMAT< '&IMAOE  FILE  NEED  TO  BE  ALLOCATED'./) 

101  FORMAT  ( 'tiREAD  THE  BLOCK  SIZE  (N  x  N)  OF  EACH  PIXEL',/) 

102  FORMAT  (ID 


